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1. Introduction 

The matter-wave diffraction has become a large field of interest over the last years, and it is extended to 
electron, neutron, atom, more massive, complex objects, like large molecules Ii-, Cqq and C70, which were 
found in experiments [1-5]. At present, There are classical and quantum methods to study interference 
and diffraction [6-12]. As is well known, the classical optics with its standard wave-theoretical methods 
and approximations, in particular those of Huygens and Kirchhoff, has been successfully applied to classical 
optics, and has yielded good agreement with many experiments. This simple wave-optical approach also gives 
a description of matter wave diffraction. However, matter-wave interference and diffraction are quantum 
phenomena, and its full description needs quantum mechanical approach. Recently, there are some quantum 
theory approach to study electron and neutron diffraction, and obtain some important and new results 
[13-17]. In viewpoint of quantum mechanics, the neutron has the wave nature, which is described by wave 
function ?/>(r, i), and the wave function ?/>(r, i) has statistical meaning, i.e., | ip(r,t) | 2 can be explained 
as particle's probability density. For the single and double-slit diffraction, if we can calculate the neutron 
wave function ip(f,t) distributing on display screen, then we can obtain the diffraction intensity, since the 
diffraction intensity is directly proportional to | ip(f,t) | 2 . In the single and double-slit diffraction, the 
neutron wave functions can be divided into three parts. The first is the incident area, and the neutron wave 
function is a plane wave. The second is the slit area, where the neutron wave function can be calculated by 
the Schrodinger wave equation. The third is the diffraction area, where the neutron wave function can be 
obtained by Kirchhoff's law. Otherwise, we consider the decoherence effect in the double slit diffraction. We 
know decoherence is the irreversible emergence of classical properties when an isolated system interacts with 
an environment [18]. The environment can be constituted by many randomly distributed particles interacting 
with the system by means of scattering processes. When these events occur in a large number, the off-diagonal 
elements of the system reduced density matrix undergo an exponential damping [19], this making the system 
to quickly lose its coherence, i.e., the decoherence is the dynamic suppression of the interference terms 
owing to the interaction between system and environment. In this paper, we study the neutron single and 
double-slit diffraction with the quantum approach, and analyze the influence of the decoherence machanism 
to the double slit diffraction. We compare our calculation results to the cold neutron (A ~ 20A) diffraction 
experiment carried out by Zeilinger et al. in 1988 [20]. We find the decoherence machanism has improved 
the calculation result of the double slit diffraction, and the theory results are agreement with the experiment 
data. 
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FIG. 1: Double-slit geometry with ai the first slit width, ai the second slit width, 6 the slit length and d the distance 
between the two slits. 



2. Quantum approach of neutron diffraction 

In an infinite plane, we consider a double-slit, its width a\ and 02, length b, thickness c and the slit-to-slit 
distance d are shown in FIG. 1. The x axis is along the slit length and the y axis is along the slit width. 
We calculate the neutron wave function in the first single slit (left) with the Schrddinger equation, and the 
neutron wave function of the second single-slit (right) can be obtained easily. At time t, we suppose that 
the incident plane wave travels along the z axis. It is 

where A is plane wave amplitude. 
The potential in the first single slit is 

fO 0<x<b,0<y< ai ,0<z<c, 
v ' y ' ' y 00 otherwise, ^ ' 

where c is the thickness of the single slit. The time-dependent and time-independent Schrddinger equations 
are 

,d ,^ , h 2 , d 2 d 2 d 2 \ 1 / —> \ 

^ + ^ + ^ + ^ ( , = o, 

where M(E) is the mass(energy) of the neutron. The relation between ip(r,t) and tl>(f) is 

4>(x, y, z, t) = 4>(x, y, z)e~i Et . (5) 
In Eq. (4), the wave function ip(x,y,z) satisfies the boundary conditions 

i/>(0,y,z) = il>(b,y,z) = 0, (6) 

ip(x, 0, z) — ip(x, 01, z) = 0. (7) 
The Eq. (4) can be solved by the method of separation of variable. By writing 



<P(x,y,z)=X(x)Y(y)Z(z). 



(8) 



The general solution of Eq. (3) is 

4>i(x,y,z,t) = ^2%l> mn {x,y,z,t) 
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Eq. (9) is the neutron wave function in the first single slit. Since the wave functions are continuous at z = 0, 
we have 

ipo(x,y,z,t) | z=0 = ipi(x,y,z,t) \ z=0 . (10) 

From Eqs. (2), (6) and (9), we can obtain the Fourier coefficient D mn by Fourier transform 

4 f 1 /*' , • n < • m7n ?^ 
Wmn = — r / / Asm — sm d^drj 

aib Jo Jo o a i 

^ m,n,odd, (n) 

otherwise, ' 

substituting Eq. (11) into Eq. (9), we can obtain the neutron wave function in the first single slit. 



Mx,y,z,t) = J2 ^mvo.mw^ sin 



16 A (2n+l)irx (2m + lWw 

sin- 



(2to+ l)(2n+ 1)tt 2 b ai 

m'n— 



2MB _ (2„ + l)2„2 _ ( 2m + l)2 x 2 ^ 

p «? e -i £ *. (12) 

The neutron wave function in the second single slit can be obtained by making the coordinate translations 
x' = x, y' = y — ai — d, z' = z, and we can obtain the neutron wave function ip2(x, y, z, t) in the second slit 

^ 16^4 . (2n+l)irx . (2m + l)-K(y-a 1 -d) 

ip 2 {x,y,z,t) = > — — — — — -J sin sin 

(2m + l)(2n + l)7r 2 b a 2 



m,n— 
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■ e V^ p ^I - V* Bt . (13) 

3. The wave function of neutron diffraction 

With Kirchhoff's law, we can calculate the neutron wave function in the diffraction area. It can be 
calculated by the formula [21] 

1 f e lkr 1 r 

ip<mt(r,t) = -— / n ■ [V'ipin + (ik - -)-ip in ]ds, (14) 

47r J s r r r 

where ip ou t(r, t) is the diffraction wave function on display screen, ipi n (r, t) is the wave function of slit surface 
(z = c) and s is the area of the aperture or slit. 
For the double-slit diffraction, Eq. (14) becomes 

1 f e ikr 1 r 

i>out{r,t) = -— / n-[V'ipi + (ik )-tpi]ds 

47r J Sl r r r 

_ 1 r £^ y/ . fc _lr 

47T J S2 r r r 

In Eq. (15), the first and second terms are corresponding to the diffraction wave functions of the first slit 
and second slit. 
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FIG. 2: Diffraction area of the single slit 

In the following, we firstly calculate the diffraction wave function of the first slit, it is 

I f e ikr 1 r 

ipou tl (r,t) = -— / n ■ [VVi + (ik - -)-ip!]ds. 

47T J Sl r r r 



(16) 



The diffraction area is shown in FIG. 2, where k 



2ME 

h 2 ' 



si is the area of the first single-slit, r' is the 
position of a point on the surface (z=c), P is an arbitrary point in the diffraction area, and n is a unit vector, 
which is normal to the surface of the slit. 



From FIG. 2, we have 



and then, 



r = R--. r '^R--r' 
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with k 2 = fcp. Substituting Eq. (17) and (18) into Eq. (16), one can obtain 

p ikR f oo oo . 

J Sn ™_n n 
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(17) 



(18) 



(19) 



Assume that the angle between k 2 and x axis (y axis) is § — a (f — and a{j3) is the angle between k 2 
and the surface of yz (xz), then we have 



k 2x — fcsina, k 2y — fcsin/3, 
n ■ k 2 = k cos 6*, 

where 9 is the angle between k 2 and z axis, and the angles 6, a, (3 satisfy the equation 



cos 2 + cos 2 (^ - a) + cos 2 (^ - (3) = 1. 



From FIG. 2, we have 



sinfi = — , 



(20) 
(21) 

(22) 

(23) 



with R = y/P + s 2 . Substituting Eqs. (20)-(23) into Eq. (19) yields 

47tk ^—^ ^— ' (2m + l)(2n + I17H 

m=0 n=0 v ' v ' 
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Jo o Jo ai 

Equation (24) is the diffraction wave function of the first slit. The neutron diffraction wave-function for the 
second slit can be obtained by making the coordinate translations x' — x,y' = y — (a + d), z' = z, it is 



4vri? ^ ^ (2m+l)(2n + l)7T 2 

m— n— 

/2M£ (2n+ 1)tt, 2 , (2to + 1)tt 
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a; dx 



" l+a2+rf e -iwV S in (2m+ 1)?r (y' - ( 0l +d))dy, (25) 

ai+d a 2 

where d is the two slit distance. The total diffraction wave function for the double-slit is 

il>out(x, y, z, t) = aipout, {x, y, z, t) + c 2 ip out2 {x, y, z, t), (26) 

where c\ and c 2 are superposition coefficients , and |ci| 2 + \c 2 \ 2 = 1. For the single-slit diffraction, we can 
obtain the relative diffraction intensity / on the display screen, 

Ioc \4> outl (x,y,z,t)\ 2 . (27) 
For the double-slit diffraction, we can obtain the relative diffraction intensity I on the display screen, 
I oc \ipout(x,y, z,t)\ 2 

= c 2 \il> outl (x, y,z,t)\ 2 + cl\tp out2 (x 7 y,z,t)\ 2 + 2c 1 c 2 Re[ip* outl (x, y, z, t)ip out2 (x, y, z,t)]. (28) 
4. Decoherence effect in double-slit diffraction 

Decoherence is introduced here using a simple phenomenological theoretical model that assumes an ex- 
ponential damping of the interferences [7, 18, 19], i.e., the decoherence is the dynamic suppression of the 
interference terms owing to the interaction between system and environment. The Eq. (26) describes the 
coherence state coherence superposition, without considering the interaction of system with external en- 
vironment. When we consider the effect of external environment, the total wave function of system and 
environment for the double-slit factorizes as [7] 

i>out(x, y, z, t) = CiVWi ® \Ei > t +c 2 ipout 2 ® \E 2 > t , (29) 

where \E\ > t and \E 2 > t describe the state of the environment. The diffraction intensity on the screen is 
now given by [7]: 

I = (1 + \a t \ 2 )(cl\ip out i{x,y 7 z 7 t)\ 2 + cl\ip out2 (x 7 y,z,t)\ 2 +2cic 2 K t Re[^* outl (x,y,z,t)ip out2 (x,y,z,t)]). (30) 
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where at =t< -E2I-E1 >t, and A t = r^fep ■ Thus, A t is defined as the quantum coherence degree. In Eq. 
(30), the two slits wave functions Vwti and ip ut2 are calculated by the quantum approach (in Eqs. (24)- 
(25)). In Refs. [7], the two slits wave functions are two Gaussian wave packets. The fringe visibility of v is 
defined as [7]: 

T T . 

(31) 



1 max ~r J-min 



+ In 

Imax and I m in being the intensities corresponding to the central maximum and the first minimum next to 
it, respectively. The value for the fringe visibility of v = 0.59 is obtained in Zeilinger et. al. experiment [20] 
(Imax — 4076, I m in = 1050), and the quantum coherence degree A t m v [7]. 

5. Numerical result 

Next, we present our numerical calculation of relative diffraction intensity. The main input parameters 
are: neutron mass M = 1.67 x 10~ 27 kg, the distance between slit and display screen I = 5m, the diffraction 
angle on yz surface a — rad, the slit thickness c — 3.0 x 10~ 5 m, the neutron energy E — 3.3 x 10~ 23 J 
(corresponding to neutron wave length A = 20A) and Planck's constant H = 1.055 x 10~ 34 Js. The equations 
(24)-(30) are series for the integer m and n. We find the series is convergence when m > 600 and n > 10, so we 
can make numerical calculation for equations (24)-(30). For single-slit experiment, the neutron wavelength 
A = 20A, the slit width ai — 90^m. In our calculation, we take the same experiment parameters above, and 
the theoretical input amplitude parameter A — 2.45 x 10 4 . From Eq. (27), we can obtain the diffraction 
intensity pattern, and it is shown in FIG. 3. In FIG. 3, the solid curve is our calculation result, and 
the dot curve is the experiment data [20]. From FIG. 3, we can find the calculation result is agreement 
with experiment data. For the double-slit diffraction, we consider two cases: coherence superposition and 
decoherence effect. For the coherence superposition, we can calculate the diffraction intensity from Eq. 
(28), and it is shown in FIG. 4. The experiment parameters are: the neutron wavelength A = 20A, the 
first and second slit width a\ = 21.9yum, a 2 = 22.5^m, the distance between the two slit d = lOOfim. In 
our calculation, we take the same experiment parameters above, and the theoretical input parameters are 
superposition coefficients c\ — 0.397, c 2 = 0.918 (|ci| 2 + | C2 1 2 = 1) and amplitude parameter A = 6.8 x 10~ 2 . 
In FIG. 4, the solid curve is our theoretical calculation, and the dot curve is the experiment data [20]. From 
the FIG. 4, we find that the theoretical result is in accordance with the experiment data, when the position 
s is in the range of |s| > 300/im. When the position s is in the range of |s| < 300/xm, the theoretical result 
has a large discrepancy with the experiment data. We find the discrepancy can be eliminated when the 
decoherence effect is considered. From Eq. (30), we can obtain the diffraction intensity pattern and it is 
shown in FIG. 5. In calculation, superposition coefficients c\ = 0.397, c 2 = 0.918, amplitude A = 6.8x 10~ 2 
and quantum coherence degree A t = 0.59. In FIG. 5, the solid curve is our theoretical calculation, and the 
dot curve is the experiment data [20]. From FIG. 5, we can find when the decoherence effect is considered, 
the calculation result is in accordance with the experiment data, and the discrepancy between the theoretical 
result and experiment data can be eliminated. 

6. Conclusion 

In conclusion, we study neutron single and double-slit diffraction with quantum theory approach. The 
calculation result of single-slit diffraction is in accordance with the experiment data. For the double-slit 
diffraction, we study the diffraction intensity by two approaches, which are the coherence superposition 
and decoherence mechanism. When we consider the coherence superposition, the theoretical result has a 
large discrepancy with the experiment data. When we consider the decoherence mechanism, the theoretical 
result is in accordance with the experiment data, and the discrepancy between the theoretical result and 
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experiment data has be eliminated. Otherwise, we think the approach has universal applicability, such as, 
it can also study electron, atom and molecular diffraction, and it can also be studied multi-slit and grating 
diffraction. 
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FIG. 3: Comparison between theoretical prediction from Eq. (27) (solid line) and experimental data taken from 

[20] (circle point) for neutron single-slit diffraction. 




(urn) 



FIG. 4: Comparison between theoretical prediction from Eq. (28) (solid line) and experimental data taken from 
[20] (circle point) for neutron double-slit diffraction, no including the decoherence effects. 
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FIG. 5: Comparison between theoretical prediction from Eq. (28) (solid line) and experimental data taken from 
[20] (circle point) for neutron double-slit diffraction, no including the decoherence effects. 



